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摘要 : 运用 REDFIT 方法 对 4 个 费 米 耀 变 体 的 MOJAVE (Monitoring Of Jets in Active 


数据 进行 了 周期 分 析 , 发 现 0851+202 


= 


Galactic Nuclei with VLBA Experiments )15. 3 GHz il 


存在 1.77 年 的 周期 成 分 1308+326 存在 8. 69 年 的 周期 成 分 1641+399 存在 2. 04 年 的 周期 


成 分 ，2251+158 存在 7. 65 年 的 周期 成 分 ， 且 这 些 周 期 成 分 的 置信 度 均 大 于 99%。 运 用 超大 


质量 双 黑 洞 模型 和 螺旋 喷 流 模型 对 这 些 周期 成 分 进行 了 解释 ， 并 估算 了 这 4 个 耀 变 体 的 中 


心 黑 洞 质量 为 1025Mo~1010Ma o 


关键 词 : 次 变 体 ，REDFIT 方法 ; 光 变 周期 ， 中 心 黑洞 质量 


中 图 分 类 号 : P157.7 ”文献 标识 码 : A 


1 引言 


HEADS fk LEA (Active Galactic Nuclei, AGN) 中 性 质 最 为 特殊 的 一 个 子 类 ， 它 


包括 蝎 虎 天 体 (BL Lac objects) 和 平 谱 射 电 类 星体 (FSRQs) 。 净 变 体 具 有 高 侦 振 和 高 光度 等 


性 质 。 炊 变 体 的 重要 特征 之 一 就 是 其 光 变 。 根据 长 期 的 观测 发 现 , 炊 变 体 不 仅 具 有 短 时 标 光 


变 , 还 具有 中 等 时 标 和 长 周期 光 变 。 利用 观测 到 的 光 变 时 标 可 以 估计 一 些 重要 的 物理 量 如 辐 


射 区 半径 大 小 、 亮 温度 、 中 心 黑洞 质量 、 辐 射 区 域 大 小 、 多 普 勒 因子 等 ， 从 而 为 辐射 模型 建 


立 提供 了 关键 特征 量 '“。 虽然 活 动 星系 核 的 光 变 曲线 一 般 是 随机 的 、 非 周期 性 的 中， 但 随 


着 近年 来 观测 数据 的 不 断 积累 和 周期 分 析 方 法 的 不 断 完善 , 越 来 越 多 的 具有 高 置信 度 的 周期 
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性 光 变 被 发 现 ， 例 如 :A0 0235-164 具有 8.13 年 的 光 变 周期 Pi，PKS 0219-164 具有 270 土 26 


天 的 光 变 周期 四 ，PKS 0405-385 具有 2.8 年 的 光 变 周期 中，S5 1044+71 具有 3. 06 士 0. 43 年 


的 光 变 周期 如，J0849+5108 具有 176 天 的 周期 中。 次 变 体 周期 性 光 变 的 研究 对 探索 活动 星系 


核 中 心 黑洞 的 结构 、 


近年 来 


物理 特性 、 动 力学 和 辐射 机 制 有 重要 作用 "™。 


国内 外 针对 活动 星系 核 周 期 性 光 变 的 研究 已 经 由 单个 源 向 数 个 源 或 大 样本 过 渡 。 


例如 ，Bhatta 等 中 对 12 颗 蝎 虎 天 体 和 8 颗 平 谱 射电 类 星体 进行 了 深入 的 时 域 分 析 ， 他 们 在 


这 些 源 的 Lomb-Scargle 和 WWZ 周期 图 中 ， 发 现 S5 0716+714、Mrk 421、PKS 2155-304, PKS 


1424-418 JU ON +325 等 具有 年 为 单位 的 准 周期 振荡 , 是 这 些 准 周期 振荡 的 置信 度 均 大 于 99%。 


Penil 等 (4 利用 费 米 大 面积 望远镜 (Fermi-LAT) 提供 的 9 年 7 射线 数据 , 使 用 十 种 不 同 的 方法 


对 2274 个 活动 星系 核 (AGN) 的 光 变 曲线 中 的 准 周期 信号 进行 了 深入 的 研究 , 发 现 具有 高 置信 


度 水 平 的 源 只 有 11 个 ， 其 中 有 9 个 源 的 周期 性 之 前 尚未 确定 ， 共 有 低 置 信 度 水 平 的 源 有 13 


个 ， 同 样 有 9 个 源 的 周期 性 之 前 尚未 确定 。 他 们 认为 活动 星系 核 周期 性 光 变 的 解释 有 多 种 ， 


例如 超大 质量 双 黑 洞 模 型 、 薄 吸 积 盘 理论 等 。 所 以 , 进一步 研究 这 类 源 的 光 变 特性 和 天 体 物 


理性 质 具 有 重要 意义 。 本 文 搜 集 了 MOJAVE (Monitoring Of Jets in Active Galactic Nuclei 


with VLBA Experiments) 的 射电 波段 观测 数据 ， 采 用 REDFIT 方法 对 这 些 炮 变 体 的 光 变 曲线 


进行 了 分 析 , 利用 


算 了 它们 的 实际 物理 驱动 周期 和 主 黑 洞 质 量 。 


超大 质量 双 黑 洞 模型 和 螺旋 吗 流 模型 对 得 到 的 光 变 周期 进行 了 解释 , 并 计 


2 观测 数据 的 选取 


MOJAVE 是 一 项 针对 北 天 区 活动 星系 核 射 电 喷 流 活 动 性 的 甚 长 基线 阵列 长 期 监测 项 目 531 


Chttps://www. cv. 


nrao. edu/MOJAVE/index. html) 。 它 的 观测 目标 有 518 个 ， 其 大 部 分 数 


据 的 观测 时 间 可 以 追溯 到 1994 年 。MOJAVE 释放 的 大 量 高 质量 数据 被 广泛 应 用 于 活动 星系 核 


的 多 波段 研究 中 。 为 了 研究 炊 变 体 可 信 的 长 周期 光 变 时 标 ， 应 选取 观测 时 间 跨 度 足 够 长 ， 数 
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据点 较 多 的 样本 。 具 体 的 筛选 条 件 如 下 


(a) 观测 时 间 跨 度 :T>15 年 ; 


(b) 年 平均 数据 点 M>4。 
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在 这 样 的 筛选 条 件 下 ， 初 步 得 选 出 了 一 批 样本 ， 但 是 为 了 得 到 的 光 变 周期 更 加 可 靠 ， 本 


文 在 周期 分 析 结 果 的 基础 上 , 给 出 了 另 一 个 限定 条 件 : 观测 时 间 跨 度 T 与 光 变 周期 P 的 比值 需 


大 于 等 于 1. 5 ， 在 此 限定 条 件 下 ， 进 一 步 简 化 样本 ， 得 到 了 满足 条 件 的 4 个 耀 变 体 源 。 表 1 


中 列 出 了 这 4 个 耀 变 体 的 基本 参数 ， 第 \1) 列 是 这 4 个 耀 变 体 的 源 名 ， 第 \2) FI Ai AAT REDE M 


的 类 型 ， 第 (3) 列 是 它们 的 Fermi-LAT 名 ， 第 \4) 列 是 他 们 的 红 移 ， 第 (5) 列 是 对 它们 观测 的 起 


止 时 间 , 98 (6) 列 是 观测 时 间 跨 度 , 第 (7) 列 是 年 平均 数据 点 数 , 用 AADP (Annual 


Points) 来 表示 。 


3 REDFIT 方法 


Average Data 


目前 ， 对 漆 变 体 周期 性 的 光 变 研究 方法 主要 有 功率 谱 方法 , Jurkevich 方法 , 小 波 分 析 


法 ， 高 散 相 关 函 数 方法 等 “1。 但 周期 提取 算法 对 时 间 序 列 的 分 析 往 往 会 受到 随机 噪声 的 


影响 ， 进 而 造成 检测 结果 的 不 确定 性 。 为 了 周期 分 析 的 准确 性 , 我 们 需要 在 周期 分 析 中 考虑 


到 随机 噪声 带 来 的 影响 以 及 对 置信 度 进 行 分 析 。 在 天 文学 、 气 象 学 和 地 理学 的 数据 中 ， 这 些 


随机 噪声 一 般 是 由 随机 过 程 产 生 的 ， 并 且 它 们 的 功率 谱 服 从 震 律 分 布 ， 


BNP (f) x 


fE Pf) 是 功率 谱 , w 是 详 指数 。 特别 地 ， 当 1 < a < 2 RUA Bri I Zr S 11 。 


红 品 声 的 存在 会 在 周期 图 的 低频 段 产生 一 系列 虚假 的 峰 ""。 


本 文采 用 一 阶 自 回 归 过 程 CART) 模型 对 燃 变 体 光 变 曲 线 中 的 红 品 声 背 景 进行 描述 站。 


ARI 模型 通常 被 用 作 零 假设 ， 以 评估 时 间 序 列 中 记录 的 变异 性 是 否 与 这 种 类 型 的 随机 起 源 


一 致 。 本 文采 用 由 Schulz 等 开发 的 REDFIT 程序 对 炮 变 体 的 光 变 曲线 进行 分 析 。 该 程 


序 直 接 从 不 均匀 间隔 的 时 间 序 列 中 估算 ART 参数 ， 即 无 需 插值 , 然后 将 估计 的 


ARI 模型 从 时 
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域 转换 为 频 域 。 其 理论 方法 是 : 对 于 任意 间隔 的 时 间 t G = 二 2,3,……N)， 有 一 个 离散 


的 ARI 过 程 r 


pi = exp(—(ti — ti-1)/7) (1) 


其 中 常数 [表示 ARI 过 程 的 特征 时 间 尺 度 ，e 表 示 均 值 为 零 ， 方 差 为 2 = 1 — exp(—2(t; — 


ti_1)/) 的 “白色 ”高 斯 噪声 。 这 个 方差 只 的 值 确保 了 AR1 过 程 的 稳定 性 并 且 使 其 具有 单位 


HŽ. AA Q) 时 域 过 程 对 应 的 频谱 Gi(fj) 为 : 


ae (2) 
Gri) = Go —) 
1 — 2p cos G +p? 


其 中 表示 截至 Nyquist 频率 yya 的 离散 频率 ，Go 表 示 平均 振幅 谱 。 平 均 自 相关 系数 p 是 由 


抽样 间隔 At = (tw — t,)/(N 一 1) 和 p = exp( 一 At/7) 的 算术 平均 值 计算 得 到 的 。 常 数 t 的 值 是 


IH Mudelsee 等 人 后 设 计 的 最 小 二 乘法 从 间隔 不 均匀 的 时 间 序 列 中 估计 出 来 的 。 在 不 需要 插 


Lu 


值 的 情况 下 ， 利 用 Lomb-Scargle 傅 里 时 变换 确定 了 一 个 不 规则 间隔 时 间 序 列 的 频谱 。 


Schulz 等 人 中 I 为 此 提出 了 一 个 计算 机 程序 ， 该 程序 利用 了 所 谓 的 


Welch-overlapped-segment-averaging (WOSA) 程 序 ， 即 将 一 个 时 间 序 列 分 割 成 重 骆 50% 的 


nso 个 片段 ， 最 终 的 光谱 估计 来 自 于 nso 个 片段 的 平均 。 有 了 Tt 的 估计 值 以 及 Go 的 适当 值 ， 就 


可 以 覆盖 公式 (2) 的 红 噪声 谱 和 从 数据 估计 的 谱 。 假 设 G- 在 每 个 频率 的 概率 分 布 服从 X 分 布 。 


还 可 以 测试 数据 频谱 是 否 与 红 噪 声 模型 一 致 。 然 而 ，Lomb-Scargle 传 里 叶 变换 在 该 方法 中 


的 应 用 是 有 缺陷 的 ， 即 与 经 典 傅 里 时 变换 相 比 ， 单 个 Lomp-Scargle 傅 里 叶 分 量 并 不 一 定 相 


互 独立 。 因 此 ， 基 于 Lomb-Scargle 变换 的 估计 频谱 可 能 是 有 偏差 的 ， 尤 其 是 ， 频 谱 高 频 端 


的 频谱 振幅 往往 被 高 估 。 因 此 ， 对 于 给 定 的 时 间 序 列 ， 基于 zt 的 无 偏 估计 的 红 品 声 谱 公式 (2) 


Tu 


不 一 定 与 同一 时 间 序列 的 “Lomb-Scargle 谱 ” 重 合 ”。 因 此 ， 此 方法 寻找 对 Lomb-Scargle 


侍 里 叶 变 换 的 偏差 校正 。 最 后 得 到 一 个 偏差 校正 的 Lomb-Scargle 周期 图 。 


chinaXiv:202211.00172v1 


ChinaXiv 合 作 期 刊 


虽然 REDFIT 程序 最 初 是 为 研究 古 气候 学 中 的 时 间 序 列 开发 的 ， 但 它 也 被 广泛 应 用 于 寻 


找 不 同 频率 下 耀 变 体 光 变 曲线 中 可 能 的 周期 性 。 例 如 Peiiil 等 人 外 利用 REDFIT 对 2274 个 活 


动 星系 核 的 光 变 曲线 进行 了 周期 分 析 ; Li 等 人 所 利用 REDFIT 对 S5 0716+714 进行 光 变 周期 


的 分 析 ， 得 到 其 约 6. 1 年 的 准 周期 ，Yang 等 人 中 | 弄 用 REDFIT 对 3FGL J0449. 4-4350 的 光 变 


曲线 进行 了 周期 分 析 ， 得 到 其 可 能 的 准 周期 振荡 为 450 天 。 


4 光 变 周期 的 分 析 


根据 MOJAVE 数据 库 中 找到 的 这 4 个 兆 变 体 源 的 射电 波段 观测 数据 ， 得 到 了 它们 的 光 变 


曲线 。 如 图 1 所 示 ( 横 轴 为 儒 略 日 ， 纵 轴 为 流量 ) 。 由 图 可 以 看 出 ， 各 个 源 的 光 变 剧烈 ,流量 


变化 出 现 了 上 升 和 下 降 的 现象 ， 且 可 能 具有 周期 性 。 


本 文采 用 上 述 的 REDFIT 程序 对 这 4 个 耀 变 体 进行 光 变 周期 分 析 。 在 REDFIT 程序 中 ,本 


文 使 用 的 输入 参数 为 n50=1( 折 县 次 数 ) ，Nsimr10000( 蒙 特 卡 洛 模拟 次 数 ) ， 


ofac=10 (Lomb-Scargle 傅 里 叶 变换 的 过 采样 因子 ) hifac-l (最 大 频率 ) 和 iwin=0( 模 型 参 


数 ) 。REDFIT 分 析 结 果 如 图 2 所 示 ， 各 子 图 中 红色 实 线 是 红 噪 声 ， 蓝 色 实 线 和 绿色 实 线 分 别 


是 置信 度 为 99% 的 X2 分布 曲线 和 置信 度 为 99% 的 蒙特 卡 洛 曲线 。 在 判断 光 变 曲线 是 否 上 共有 周 


期 信号 时 ， 本 文 根 据 光 变 曲 线 的 周期 图 中 的 峰值 是 否 同时 达到 99% 的 x* 和 99% 蒙 特 卡 洛 置 信 


度 水 平 ， 并 结合 光 变 曲线 的 变化 特征 和 红 噪 声带 来 的 影响 (在 周期 图 的 低频 段 产 生 一 系列 虚 


假 的 峰 ) 进行 判断 。 周 期 图 中 周期 信号 对 应 的 频谱 峰值 满足 这 个 条 件 时 ， 我 们 认为 它 


能 的 周期 信号 ,如 果 这 个 周期 信号 对 应 的 振荡 次 数 在 3 次 以 上 , 我 们 则 认为 这 是 一 个 强 周 


期 信号 。 在 周期 图 中 ,特别 是 高 频 部 分 可 能 会 存在 部 分 伪 周 期 信号 ， 这 可 能 是 由 于 观测 数据 


的 采样 率 较 低 ，REDFIT 程序 不 能 分 辨 造成 的 。 
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2(A) 是 源 0851+202 的 周期 图 ， 从 图 中 发 现 它 具有 约 1.77 年 的 强 周期 信号 。Yuan 等 


PVE TE 12. 02 土 0. 41 年 的 光 变 周期 。 图 2 (B) 是 源 1308+326 的 周期 图 ， 从 图 中 可 以 发 现 


它 具 有 8. 69 年 的 强 周 期 信号 。Qian 等 所 获得 的 16. 9 土 0. 3 的 光 变 周期 可 能 是 这 个 8. 69 年 


的 周期 成 分 和 加 2 次 而 成 的 。 图 2 (C) 是 源 1641+399 WAWE, 从 图 中 可 以 发 现 它 具 有 2. 04 


年 的 强 周期 信号 。Berton 等 "9 获得 了 它 0. 28 年 的 光 变 周期 。 图 2(D) 是 源 2251+158 的 周期 


图 ， 从 图 中 可 以 发 现 它 具 有 7. 65 年 的 强 周 期 信号 。Das 等 中 "获得 了 9 年 的 光 变 周期 ， 由 于 


样本 之 间 存 在 差异 ， 两 个 周期 值 的 关系 难以 判断 。 表 1 第 (8) 列 给 出 了 我 们 得 到 的 这 些 强 周 


期 信号 。 由 此 可 见 ， 粮 变 体 的 周期 性 光 变 是 十 分 复杂 的 ， 周 期 成 分 可 能 并 不 唯一 。 
表 1 READ AMES AS 


Tab. 1 Possible periodicities in the radio | ight curves and central black hole masses 


of blazars 
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Source Source Fermi-LA Redshi Observati Time AADP Period yp Pa M 
Name Class T Name ft on Time Span 
(yr) (yr) (yr) (Mo) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
08514202  BLLac J0854.8+ 0.30 1995.4.7-20 26.7 6. 65 1.77 11.84 188.71 5.6X10°+1 
2006 6 21.12.10 1.7 士 3.92  .1x10' 
1308+326 FSRQ J1310.5+ 0.99 1995.1.20-2 26.67 4.04 8. 69 14.64 927.43 7.1X10°+ 
3221 73 021.9. 13 0.7 士 2.13  4.3x10 
1641+399 FSRQ J1642.9+ 0.59 1995. 1. 20-2 26.67 4. 46 2. 04 20.8+ 554.04 3.1X10" 土 
3948 3 021. 9. 13 1.0 士 1.28  1.9x10 
2251+158 FSRQ J2253.9+ 0.85 1995. 5.27-2 26.56 4.65 7. 65 15.74 1014. 3 8. 2x10" 
1609 9 021. 12. 10 1.0 324.1 a 2X 
2 10’ 


- 
76 ne^ Ma 
~ 


图 1 MOJAVE 15.3 GHz 光 变 曲线 


Fig. 1 Radio light curves of blazars at 15.3 GHz obtained from the MOJAVE program 
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线 进行 分 析 


Fig.2 The resulting periodicity analysis of the radio light curves of blazars using the REDFIT method 


5 讨论 和 结论 


本 文 收集 MOJAVE 数据 库 中 


的 射电 波段 观测 数据 ， 利 用 


REDFIT 对 0851+202 等 4 个 耀 变 


体 源 进 行 了 光 变 周期 的 分 析 , 在 0851+202 的 光 变 曲线 中 发 现 它 存在 1. 77 年 的 周期 成 分 , 在 


1308+326 的 光 变 曲线 中 发 现 它 存在 8. 69 年 的 周期 成 分 ， 在 1641+399 的 光 变 曲线 中 发 现 它 


存在 2. 04 年 的 周期 成 分 ， 在 2251+158 的 光 变 曲线 中 发 现 


吸 积 盘 的 细 盘 模型 59。1980 年 ，Begelman 等 中 第 一 次 


它 存 在 7. 65 年 的 周期 成 分 。 


目前 ， 周 期 性 光 变 的 理论 解释 主要 有 :超大 质量 双 黑洞 模型 ”下 、 螺 旋 喷 流 模型 中 1、 


是 出 了 超大 质量 双 黑 洞 的 概念 ， 认 为 


双 黑 洞 是 星系 合并 产生 的 。 双 黑洞 模型 多 用 来 解释 光 变 周期 的 起 源 ， 该 模型 认为 黑洞 的 轨道 


运动 引起 了 吸 积 盘 和 喷 流 的 进 动 , 喷 流 和 视线 方向 的 夹 角 变 化 引起 了 观测 上 的 光 变 。 周 期 性 


光 变 、 扭 曲 的 喷 流 和 发 射线 的 双 峰 结构 成 为 超大 质量 


等 BI 和 Valtonen 等 5 用 超大 质量 双生 


等 BI 和 Li 等 69 用 超大 质量 双 肌 


EE 双 黑 洞 系统 存在 的 间接 观测 证 据 中 1。 

许多 研究 都 使 用 超大 质量 双 黑 洞 模 型 对 长 时 标 光 变 的 周期 性 进行 了 解释 ， 例 如 :Sillanp 音 
黑洞 模型 解释 了 0J287 约 12 年 的 准 周期 振荡 。Tavani 
a 洞 模型 解释 了 其 他 源 的 光 变 周期 。 陈 军 平等 用 超大 质量 双 


黑洞 模型 解释 了 耀 变 体 CGRaBS J0929 +5013 和 J2146-1525 KJA. AZRA 


也 用 超大 质量 双 黑 洞 模型 解释 了 耀 变 体 CGRaBS J0835+6835 的 光 变 周期 。 
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本 文 使 用 超大 质量 双 黑 洞 模型 和 螺旋 喷 流 模型 对 粮 变 体 光 变 周期 进行 解释 。 在 超大 质量 


双 黑 洞 模型 中 ， 观 测 准 周期 P 与 实际 物理 驱动 周期 Pr 之 间 的 关系 为 上 


Yb (3) 
Py = P 3 
4 qum 


其 中 红 移 z 在 表 1 中 已 经 给 出 ;加 为 体 洛 伦 效 因子 ， 各 耀 变 体 源 的 体 洛 伦 效 因子 已 经 列 入 表 


1 第 (9) 列 中 中， 将 使 用 REDFIT 得 到 的 周期 代入 (3) 式 得 它们 各 自 的 实际 物理 驱动 周期 Pa， 


如 表 1 第 (10) 列 所 示 。 如 果 主 黑洞 和 次 黑洞 之 间 的 质量 比 R <<, 称 之 为 “ 主 并 合 ”; 如 果 


3 < R<104， 称 之 为 “次 并 合 ””。 不管 是 哪 种 并 合 ， 主 黑洞 的 质量 计算 都 可 以 用 下 列 公 


itti $m 


3 
M ~ P? R5105Mg (4) 


其 中 ， 对 于 超大 质量 双 黑洞 系统 的 主 次 黑洞 质量 ,假设 它 为 R = 2 07. ER 以 及 各 个 次 变 体 


源 的 实际 物理 驱动 周期 代入 (4) 式 得 到 它们 各 自 的 主 黑 洞 质量 ， 结 果 见 表 1 8 0D Fi]. EAE 


体 的 光 变 由 喷 流 主导 , 观测 到 的 周期 可 能 起 源 于 喷 流 的 螺旋 进 动 。 而 喷 流 的 螺旋 进 动 由 超大 


质量 双 黑 洞 的 轨道 运动 驱动 , 也 就 是 说 用 超大 质量 双 黑 洞 模型 来 解释 这 些 炮 变 体 的 中 长 时 标 


光 变 是 可 信和 的 。 


此 外 ， 螺 旋 结 构 的 喷 流 模型 也 会 导致 Blazar 天 体 的 周期 性 光 变 。 由 强 磁场 或 轨道 运动 


导致 的 扭曲 喷 流 ， 在 视角 改变 时 ， 由 于 多 普 勤 效应 会 使 得 观测 的 流量 产生 变化 。 对 于 非 均匀 


旋转 螺旋 喷 流 ， 其 中 发 射 物质 速度 的 方向 与 我 们 的 视线 夹 角 会 周期 性 的 变化 。 由 于 多 普 勒 增 


亮 效 应 ， 就 会 导致 观测 的 流量 周期 性 变化 中。 这 4 个 源 都 是 炮 变 体 ， 其 辐射 流量 主要 来 自 


于 喷 流 中 的 非 热 辐射 ， 而 超 长 基线 阵列 (VLBA) 观测 发 现 其 喷 流 呈现 扭曲 的 形状 中， 所 以 其 


周期 性 光 变 也 可 能 是 螺旋 喷 流 旋转 造成 的 。 
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Quasi-periodic Variations in the Radio 
Light Curves of Fermi Blazars 


Wen Yi, Xiao Yuting, Li Xiaopan, Yang Haiyan 


(School of Physics and Information Engineering，Zhaotong University，Zhaotong 
657000) 

Abstract: We search for possible periodicity in the MOJAVE (Monitoring Of Jets 
in Active Galactic Nuclei with VLBA Experiments) observations at 15.3 GHz of eight 
Fermi blazars using the REDFIT method. Possible periodic components of about 1. 77, 
8.69, 2.04 and 7.65 years were found in the sources 0851+202, 1308+326, 1641+399 
and 2251+158 at a 99% confidence level, respectively. The super-massive binary black 
hole model and the helical jet model have been used to explain these 
periodicities. And they are also used to caculate the central black hole masses of the f 
our sources which is 10?M-10!? Mo. 

Key word: Blazar; REDFIT method; Period of light variation; Central black hole 


masses 


